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Abstract 

We generalize modern ideas about the duality between Wilson loops and scat- 
tering amplitudes in J\f = 4 SYM to large N QCD by deriving a general relation 
between QCD meson scattering amplitudes and Wilson loops. We then investigate 
properties of the open-string disk amplitude integrated over reparametrizations. 
When the Wilson loop is approximated by the area behavior, we find that the QCD 
scattering amplitude is a convolution of the standard Koba-Nielsen integrand and a 
kernel. As usual poles originate from the first factor, whereas no (momentum depen- 
dent) poles can arise from the kernel. We show that the kernel becomes a constant 
when the number of external particles becomes large. The usual Veneziano ampli- 
tude then emerges in the kinematical regime where the Wilson loop can be reliably 
approximated by the area behavior. In this case we obtain a direct duality between 
Wilson loops and scattering amplitudes when spatial variables and momenta are 
interchanged, in analogy with the AA=4 SYM case. 



*Also at the Institute for Advanced Cycling, Blegdamsvej 19, 2100 Copenhagen 0, Denmark 



1 Introduction 



The relation between planar diagrams and dual resonance models has a long history since 
the pioneering works [1]. A long-standing belief [2] is that SU(N) Yang-Mills theory is 
equivalent at large N to a free string, while the 1 /iV-expansion corresponds to interactions 
of the string. 1 A great recent progress along this line is associated for J\f = 4 super Yang- 
Mills theory (SYM) with the AdS/CFT correspondence [4, 5] (see Ref. [6] for a review), 
where the strong-coupling limit of SYM is described by supergravity in anti-de Sitter 
space AdS$ x S* 5 . 

The (finite part of the) 4-gluon Maximally-Helicity- Violating on-shell scattering am- 
plitude in SYM theory has the form 

A(s,t) = A tIcc e^ l °^ (1.1) 

(where s and t are usual Mandelstam's variables) as was conjectured [7] on the basis 
of three-loop calculations. To explain Eq. (1.1), the Wilson-loop/scattering-amplitude 
(WL/SA) duality was introduced [8] at large 't Hooft couplings A, which has been then 
advocated in SYM perturbation theory [9]. This duality (for a review see Ref. [10]) states 
that the scattering amplitude (divided by the kinematical factor A tree ) equals the Wilson 
loop for a polygon whose vertices X{ are related to the momenta pi of scattering gluons 
by 

Pi = K (xi - , (1.2) 

where K = Xjl-noi is the string tension. 

Our goal in this paper is to find out what features of the described WL/SA duality 
(if any) remain valid for QCD and, in particular, how is it possible to maintain the 
relation of the type (1.2) which would relate large momenta in scattering amplitudes 
with loops of large size. Of course this is not possible in QCD perturbation theory, 
where \p\ ~ l/|x| because of dimensional ground. But nonperturbatively a dimensional 
parameter K ps (400 MeV) 2 appears in QCD, which shows up in the area-law behavior 
of asymptotically large Wilson loops: 

W(C) lar S C e-**-<°>, (1.3) 

where S m i n (C) is the area of the minimal surface bounded by C, that results in confine- 
ment. Strictly speaking, this requires the limit of the large number of colors N or the 
quenched approximation. 

As is well-known by now, a string theory, which QCD is supposedly equivalent to, 
is not the simplest Nambu-Goto string. Some extra degrees of freedom living on the 
string are required which are most probably conveniently described by a presence of extra 
dimensions. The asymptotic behavior (1.3) is nevertheless universal for large loops. 

In this paper the main topic is an investigation of the relation between the Wilson 
loop and the corresponding amplitude in large N QCD, or, alternatively, quenched QCD 

1 See e.g. Ref. [3] for an introduction and review of the old works on the QCD/string correspondence. 
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for any N. In general this involves integration over an infinite number of loops. However, 
if one considers large loops there is a considerable amount of evidence from lattice gauge 
calculations in 2+1 and 3+1 dimensions for various iV's that the Nambu-Goto action 
describes the behavior of the Wilson loops quite well. To give an example, in [11] it is 
shown with unprecedented precision that the static quark potential in quenched SU (3) 
lattice gauge theory is well described by the first two terms (oc r and the Liischer term 
oc 1/r) in a long-distance expansion of the Nambu-Goto action. Furthermore, in both 
three and four dimensions the transition from perturbative to string behavior takes place 
"at surprisingly small distances" [11]. There also exists a number of other comparisons 
between results from the Nambu-Goto action, e.g. the closed string spectrum and SU (N) 
for various values of N, see [12] where further references can be found. 

The various results from lattice gauge theories can be summarized by the statement 
that the Nambu-Goto action describes the large and not so large [11] distance behavior 
of quenched QCD surprisingly well. This action has the well known anomaly for d ^ 26, 
which however is suppressed for long strings [13]. The remarkable success of the Nambu- 
Goto string as an effective action led us to reconsider the relation between the Wilson 
loop W(C) and the corresponding amplitude. The idea is to perform the unpleasant sum 
over all C's by inserting the Nambu-Goto action in W(C). Then the sum over C becomes 
an integration over the string field x^(cr, r), and is at least in principle controllable. In 
practice, in this paper we start out with a more modest program, where only the area 
behavior of W(C) is inserted in the form of the disk amplitude. Thus, the effects of the 
Liischer term is not included in the investigation reported in this paper. 

The present paper is an extended version of Ref. [14], but some of the reported results 
are novel. In Sect. 2 we derive a relation between the M-meson scattering amplitude and 
the Wilson loop for the case of fermion and scalar quarks, valid for QCD in the large 
N limit. Our results are general in the sense that if the Wilson loop is known one can 
obtain the scattering amplitudes by performing some path integrations. In Sect. 3 we 
then take up the old idea that the large N QCD Wilson loop should be identified with 
the disk amplitude in certain string models as far as the leading large distance behavior 
is concerned. We emphasize that it is crucial in this construction to integrate the string 
disk amplitude over reparametrizations of the boundary contour. In Sect. 4 we show, 
taking a functional Fourier transform to momentum space, how the disk amplitude leads 
to the Koba-Nielsen amplitude. Integrating over the reparametrizations, we also derive 
project ive-invariant off-shell scattering amplitudes. 

In Sect. 5 we return to the general formula from Sect. 2 for the relation between the 
meson (made from fermion quarks) scattering amplitude and the Wilson loop. Here we 
insert the area law (with no sub leading perimeter term, Liischer term, ...) and for the 
(off-shell) M-particle amplitude we derive the formula 



M-l ^i 1 m r 
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(1.4) 



where the Ap's are particle momenta and where the "kernel" /C in general has no mo- 
mentum dependent singularities for fa — > fa. The first factor in the integrand on the 
right-hand side is the well-known Koba-Nielsen integrand, and it produces poles in the 



representation of the Wilson loop, which means that the area should be large and cor- 
respondingly the momenta should also be large. Thus, due to this condition there is no 
tachyon (or other low lying states) in the spectrum, as one would indeed expect in QCD. 
We further show that in the case of large M the kernel K, degenerates to a factor which 
is essentially independent of the Ap's. This implies that when the area behavior of the 
Wilson loop dominates the dynamics we get the interesting result that when many par- 
ticles are produced in some collision then the scattering is given by a Veneziano type of 
amplitude. Of course, this is for large N QCD, but it would be interesting to see to which 
extent this would be valid at the LHC collider with N = 3. 

The condition that the various momentum transfers should be large in order that the 
area behavior for the Wilson loop can be inserted is only a necessary but not sufficient 
condition for the validity of our approach. This is important in view of one of the histor- 
ical reasons for not using the Veneziano amplitude in strong interaction phenomenology: 
at large transverse momenta this amplitude decreases exponentially in contrast to the 
experimental data. If we consider the 4-point function in terms of the usual Mandelstam 
variables s and t, our approach is valid when 1/a' < — t <C s in Minkowski space, be- 
cause it dominates over other contributions. However, if —t ~ s, this is no longer true, 
because then the Veneziano amplitude becomes a tiny exponentially decreasing function 
and other contributions lead to a power decrease of the amplitude, which of course is far 
more important than an exponential decrease. Therefore the area behaved Wilson loop 
no longer dominates when —t ~ s. 

In Sect. 6 we find a relation between the space variables and the external momenta 
which is the QCD analogue of the SYM WL/SA duality exhibited in Eq. (1.2). If z{fa) 
describes the contour of the Wilson loop as a function of an angular parametrization 0, 
then 



where pi are related to the external momenta Api by Api = — pi, and is the usual 
step function. If this expression is inserted in the z-dependent Wilson loop, it reproduces 
the scattering amplitude with a large number M of external particles when the fa's are 
integrated over. Thus, the relevant contours are given by the constant momentum vectors, 
in complete analogy with Eq. (1.2). The difference with the supersymmetric case is that in 
QCD we have to integrate over the parameters fa which represent the points in parameter 



integral for fa 



fa. Physically these poles only occur when the area law is a good 




(1.5) 



i=i 
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space where the momenta enter. Otherwise (1.5) is like Eq. (1.2), since z(4>) equals Pi/K 
in the interval from the point X{ to Xj+i, so the vector Xj+i — X{ equals Api/K, as in 
Eq. (1.2). 

Some more technical details are discussed in the Appendices. In Appendix A we give 
an example of how to operate with path integrals related to the ordering of gamma ma- 
trices. In Appendix B we review the modern approach to the minimal area as a boundary 
functional. To illustrate the asymptotic area behavior, we evaluate in Appendix C the 
integral over reparametrizations in the disk amplitude for a large circle. Appendix D is 
devoted to calculations of the path integral over reparametrizations. In Appendix E we 
consider a scattering amplitude which appears in QCD, when the Wilson loop is sub- 
stituted by an exact area law, i.e. Eq. (1.3) holds not only asymptotically but for all 
contours. 



2 QCD amplitudes dual to Wilson loops: general re- 
sults 

In large N QCD, Green's functions of M colorless composite quark operators (e.g. q(xi)q(xi)) 
are given by the sum over all Wilson loops passing via the points Xi (i — 1, . . . , M), where 
the operators are inserted. This approach was first advocated on the lattice [15] and then 
extended [16] to the continuum. To obtain a scattering amplitude, one makes the Fourier 
transformation with respect to Xj and takes the corresponding momenta Pi on shell. 

The weight for the summation over paths depends on both the quark spin and the 
quantum numbers of the operators (see Refs. [17, 3] for more detail). The simplest results 
are for fermion quarks and the scalar operators q(xi)q{xi) when the connected Green's 
function is 

G(x u ...,x M ) = (jJq(x i )q(x i )\ = £ W(C) . (2.1) 

\i=l I conn C3xi,...,xm 

Here W(C) is the Wilson loop in pure Yang-Mills theory at large N. For finite N, 
correlators of several Wilson loops have to be taken into account. 

In Euclidean space the standard weight for the summation over paths in Eq. (2.1) 
reads explicitly 

m \ m °2 

* = 1 / conn * = 1 o 

V Zi (t) J Vk(t) sp Pexp j i J dt[ Zi (t)k(t) - ~f(t)k(t)\ J W(C), (2.2) 



z i (0)=x i _ 1 



where the i-segment of the loop C from x^i to X; L is represented by the function Zi(t) 
(0 < t < T,j) and x M = (since the loops are closed). The integration over the variable 
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gamma matrix 7^(t) is explained in Appendix 7. Above we have introduced the convention 
to be used in the rest of the paper that dot means derivative with respect to whatever 
argument a function has. Thus, for some function / we have f(x) = df(x)/dx, whereas 
f(y) = df(y)/dy. 

Equation (2.2) is essentially derived in [3], where further references can be found. For 
the readers convenience we repeat the essential steps. In QCD the quark fields can be 
integrated out and we have the Feynman disentangling 2 



1 

\x 



oo 

J dr (y\e- T ^» +m) \x) 



liNn + m 

u 

oo 

= J dr e~ mT J Vk„ sp P e -/o T dt7M(t)(iM*)-i^M) s(k fl (t) + id^t)) S(x - y), (2.3) 
o 

where sp and P act on the gammas and color matrices. The functional delta function has 
the representation 

Silent) + id„(f)) S(x ~y) = J Vv » e^MfXMO+iM*)) $( x - y). (2.4) 
Here f M (t) has no restrictions. Then 

5(k^t) + id^t)) S(x - y) = J Vv^e^o dtv ^ k ^ S(x + dtv(t)-y). (2.5) 
With Zn(t) = x >1 + f*dt'v fl (t') we get 

6(k li (t)+id li (t))6(x-y)= J Vz li e i K dti »® k »®. (2.6) 



z(0)=x 
z(r)=y 



Introducing Eq. (2.6) in Eq. (2.3), we obtain 



oo 

/ dT< 



y\e- T ^» +m) \x) 



oo 



J dr e - mT J Vk„ J Vz„ sp P e 1 ^ dt l^-j^+j^] 



(2.7) 



z(0)=x 
z{r)=y 



Shifting the integration by + A^, we get 

dr(y|e- r( > v " +m) |:r) 



oo 

= J dr e~ mT j Vk^ J ^spPe 1 ^*^-^^]. ( 2 .8) 



z(0)=x 
z{r)=y 



2 We use the notation of [3]. The states \x) etc. are eigenstates of x and the standard Feynman 
disentangling is used in (2.3). 
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This result leads immediately to Eq. (2.2) since the contour of C is composed of pieces 
each of which can be represented as in Eq. (2.8). 
Introducing new proper-time variables 

i 

% = Y^r 3 , T = T M (2.9) 

we rewrite Eq. (2.2) in a time-ordered form: 
G (xi, . . .,x M ) 

°° M-l r f 

= / dTe- mT J] d% / Vz(t)Vk(t) sp P e^^WW-TWfcW] ^[^)] , 

4=1 *(o)=x 

z ( r i)= a: i 

z(^)=^A'/= a: 

(2.10) 

where the loop is represented by a single function z^(t) (0 < t < T): 

C 3 X!,...,x M = {z^(t) , z(0) =x ,.. .,z(7l) = Xi, . . .,z(T) = x M = x } , (2.11) 
describing a closed loop passing via the (ordered) set of points X{. It is essential here that 

Finally, we introduce the angular-type variables 

= 2tt^ (0 < < 2tt) , 

<t>i = 2vr^ (O<0,<0 W <2vr), 
0m = 2tt (2.12) 

so that 



oo 

G(xi, . . .,x M ) — J dT (i^j 



M-l 

„ -mT 



M-l 

X 

i=l 







JI / y 2M0) y m(0) SP p js?wmm-T«m*)/**\ w[z{(f))] . 



z(0)=z(27r)=x 
z (<t>i)= x i 



(2.13) 



The on-shell M-particle amplitude can be obtained from the Green function (2.13) by 
applying the standard LSZ reduction formula. When making the Fourier transformation, 
it is convenient to represent M momenta of the (all incoming) particles by the differences 

A Pi = Pi _ 1 - Pi . (2.14) 
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Then momentum conservation is automatic while an (infinite) volume V is produced, say, 
by integration over x$. We therefore define 

i M r 

G (A Pl , A PM ) = 77 II / d ^ eiE< APiXiG {x!,...,x M ). (2.15) 

i=i "* 

We can further rewrite this formula, introducing a momentum-space loop p^((f>) which 
is piecewise constant: 

p(4>) = Pi for <f>i < <f) < <fi m . (2.16) 

Noting that 

m = -J2 ^Pi K<t> - 00 , Apl = pi_, - pi , (2.17) 

i 

we write 

^ Ap^ = - / d0p(0) • z(0) = / d0p(0) • i(0) (2.18) 
j ./o ./o 

which is manifestly parametric- invariant. 

Inserting (2.13) into Eq. (2.15) and noting that 

n / ^ / ^(^) p ^)) = J vz ^) p ^)) > ( 2 - 19 ) 



i=i 



z(0)=z z(0)=z(2tt) 

z(4> i )=x i 



where F is some functional, we obtain 

<f>m 

G(A Pl ,...,Ap M ) = / dT 



2?r e" mT n / d ^ / ^(0) 

x / vz((j>) sp p ^C^w+pmm-rim^)/^] W [z{<t>)\ , (2.20) 

z(0)=z(2vr)=0 

where p(<f>) is piecewise constant as is given by Eq. (2.16). We do not integrate over 
z(Q) = z{2n) which would produce the (infinite) volume factor because of translational 
invariance. 

In the case of scalar quarks the above procedure can be repeated and we obtain for 
the amplitude in position space 

G(x 1 ,...,x M ) = (T[<pt(x i )<p(x i )\ = W ^ C ) ( 2 - 21 ) 

\i=l / conn CBxi,...,x M 

or, more explicitly, 

M-1 



M-1 9,1 |J 

X 



fljtyi J VzWe-S^WMWiztt)], (2.22) 

i = l n ~^n\ ./o~\ 



3(0) = Z (27r)=a: 
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where /i — 2ir/T. In momentum space this gives 

1 °f /r\ M " 1 
G(Ap u ...,A PM ) = ±JdT^ e-"^ 

o 

M-l ^ 

x[]|d^ | ^(0)e-^" d ^ i2 W +i ^ d ^)-W^[ 2 (0)]. (2.23) 

1=1 O z(0)=z(2tt)=0 

In this case there is no A;— integration and no term involving the gamma matrices. The 
z 2 — term is specific for the scalar case. It is related to the occurrence of the second 
derivative in the Klein-Gordon operator. Introducing an auxiliary k— integration, we can 
rewrite Eq. (2.23) in the form 

1 r / t \ r r 

G (A Pl , A PM ) = — J dT y^J e U J d ^ J Vk ^ 

o 1=1 o 

Vz{<t>) e i / ( r^[(fcW+fW)-W-^ 2 W/4-]^[ z (0)] (2>24) 

2(0)=2(2tt)=0 

which looks similar to Eq. (2.20). 




3 Wilson loop as string disk amplitude 

An old idea is to identify the Wilson loop of large iV QCD with the (tree level) disk 
amplitude in a certain string theory which QCD is equivalent to. As is already pointed 
out, the simplest Nambu-Goto string in flat space appears to be surprisingly accurate for 
large loops and reproduces the asymptotic area law (1.3). We shall therefore first review 
the known results for the bosonic string. 

The calculation of the tree level disk amplitude for the Polyakov string has a subtlety 
associated with fixing conformal gauge [18, 19, 20, 21]. The decoupling of the Liouville 
field <f(r,a) is possible only in the interior of the disk, while its boundary value <p(l,a) 
determines the metric at the one-dimensional boundary: 

h(a) = e^ a)/2 . (3.1) 

At the classical level this fixes the parametrization of the boundary contour. The path 
integral over the boundary value of the Liouville field then restores an invariance under 
reparametrizations of the boundary in quantum case. 



3.1 Unit circle parametrization 

Let us parametrize the unit disk D by the variables < r < 1 and a G [0, 2tt). Alterna- 
tively, one can map the unit disk onto the upper half-plane: 

1 + r e iCT . . 

z = i-— -. 3.2 

l-re la v ' 
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Then the boundary <9D, associated with r = 1, is mapped onto the real axis parametrized 
by — oo < s < +00, so that 

s(a) = — cot — . (3.3) 
2 

Integrating over the string fluctuations inside the disk, i.e. over X(r,a) with r < 1, 
we arrive at the boundary action which is a functional of the field 

X(l,a)=x(a) (3.4) 

at the boundary. 

Actually the original Fradkin-Tseytlin calculation [22] used the disk parametrization. 
Their result for the disk amplitude of bosonic string in d = 26 implies 

V[x(a)} = exp ( - j J d<7id<7 2 x (a 1 )G~ 1 (a, - a 2 ) x(a 2 )^j , (3.5) 



where 



K = (3.6) 



is the string tension and 

1 1 00 1 

G(a) = — — log (2 — 2 cos a) — — — cosmcr (3.7) 

m=l 

is the corresponding Neumann function for the disk. It obeys 

d 2 1 1 °° 

G^ 1 (a) = — -— -G(cr) = ; - = — > m cos ma . (3.8) 

y 1 da 2 y ' 2n(l -coscr) vr ^ y ' 

v ' m=l 

Using Eq. (3.8) we can rewrite the disk amplitude (3.5) as 

*[f (<t)] = exp ^ - y J daxdfTs f ((Ti)G (<Ti - <7 2 ) f((7 2 )^ . (3.9) 



This formula is pure classical: the exponent is just the classical boundary action while 
determinants coming from the integration over the fields inside the disk are ignored. For 
this reason the same result holds for the bosonic part of the classical boundary action of 
superstring in d = 10. But more subtle phenomena, such as the Luscher term, which are 
due to determinants are not captured by Eq. (3.9). 

3.2 Integration over reparametrizations 

The exponent in Eq. (3.9) is not invariant under reparametrization of the contour: 

x(a) ->£ (0(a)), (3.10) 
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where the reparametrizing function 9(a) obeys 

9(0) = 0, 9(2n) = 27r, ^1>0. (3.11) 

da 

Polyakov [23] proposed (see also Ref. [24]) to integrate the disk amplitude (3.9) over 
the reparametrizations thus providing a reparametrization-invariant disk amplitude which 
can be identified with the Wilson loop in large N QCD. For a unit-circle parametrization 
this gives 

*[£(•)]= J V diS 9(a)exp(^-^ J da,da 2 x (9(a,)) G' 1 (a, - a 2 ) x (9(a 2 )) )j 

o 

(3.12) 



Diff(S^) 

or, equivalently, 



2tt 

*[£(•)]= J V diS a(9)ex P ^-jJd9 1 d9 2 x*(9 1 )G(a(9 1 )-a(9 2 ))x*(9 2 )^ . (3.13) 

DiffiS 1 ) 

Here the path integration is over the (infinite) group of diffeomorphisms of a circle 
Dif f(S l ), i.e. over functions with non-negative derivative 9'(a), which obey Eq. (3.11). 

An explicit expression for the measure for integrating over reparametrizations can be 
given using an expansion over the complete set of basis functions fj(a): 

L 

W = (3-14) 

where fj = 1 at the j-th interval [(Jj-i, &j] (er = a L — 2n) with arbitrary a^s and vanishes 
otherwise. A continuous function is approached when the number L of intervals becomes 
infinite and aj — aj-i — > 0. We can then define the measure V d is 9(a) by 

I v„ m .. . = | n / • • • ■ (3.i5) 

Diff(Si) ^ 

If a is chosen to be the length of an arc of the unit circle, then aj = 2nj/L. Otherwise, we 
have <7j = a (2itj/L). As distinct from the usual measure T>9(a): the integrals in Eq. (3.15) 
are ordered, while the additional factors are needed to provide necessary symmetries as 
is explained in Sect. 4. 

If the unit circle is mapped onto the real axis by Eq. (3.3), the disk amplitude takes 
the form 

+oo 

*[*(•)]= / Pdtt*(g)exp ( - ^ J ( ^ d * 2 2 )2 [x(t(si)) -^(* 2 ))] 2 ) , (3.16) 

Diff(R) -oo 
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where the path integral over t(s) restores the invariance under reparametrizations 

s^t(s), i(-oo) = -oo, t(+oo) = +oo, ^>0. (3.17) 

as 

The measure on Diff(M) can be given by 

+00 L _ 1 tj+i 

J ^ i(s )....n,J d(i |^n/^|^---, (3.18) 

Diff(R) -00 3 -"--oo 

which is of the same type as in Eq. (3.15). 

The equivalence of Eq. (3.16) and Eq. (3.12) can be shown using Eq. (3.3). For the 
real-axis parametrization the Green function and the inverse one read 

G(s) = ——log |s| (3.19) 

7T 

and 

G-\s) = -—. (3.20) 
%s z 

These formulas are the counterparts of Eqs. (3.7) and (3.8), while that of Eq. (3.13) is 



+00 

*[£(•)]= J Pdiffs(t)exp^ y d*idf 2 f(fi) log 1 8(h) -s(t 2 ) |f(* 2 )^ . (3.21) 

Diff(R) -oo 

To show the equivalence of (3.16) and (3.21) (or (3.12) and (3.13)), we integrate the 
exponent by parts 

+00 +00 



— oo 

+00 

= \ J d tl dt 2 [X(h) - X(h)} 2 ^g [8(h) ~ S(t 2 )] 2 

— oo 
+00 

= - J d*idt 2 £(ti)£(* 2 )log[s(*i)-s(f 2 )] 2 . (3.22) 



3.3 Large loops and minimal area 

In spite of the fact that the right-hand side of Eq. (3.16) is derivable, as is already 
mentioned, for bosonic string in d = 26 or superstring in d — 10, where the integration over 
reparametrizations comes from the integration over the boundary value of the Liouville 
field, we shall use it only as an ansatz for asymptotically large loops or, equivalently, very 
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large K, when the integral over reparametrizations has a saddle point at t(s) = t*(s). 
This will be crucial for reproducing the area-law behavior 

*[£(.)] lars «° ops e~ KS ™^ (3.23) 

for asymptotically large loops. 

The appearance of the area of the minimal surface, spanned by the contour 
minimum of the boundary action in Eq. (3.16) (or Eq. (3.12)) is related to the fact that 
it is nothing but the functional known in mathematics as the Douglas integral [25] , whose 
minimum with respect to reparametrizations does give the minimal area: 

,[#))-«)f 



min '« \^ J ds J ds 



(s - s'y 



+oo +oo 



- j ds j ds' L V V (;_ S ; )2 V = ^ [*(•)] • (3.24) 



— oo — oo 



This issue is clarified in Appendix B. 

The necessity of the function t*(s) (or 9*(cr)), reparametrizing the boundary is due 
to the fact that coordinates describing the minimal surface have to obey the condition of 
the conformal gauge for the quadratic action, used in the Polyakov string formulation, to 
coincide with the area. In particular, Eq. (3.4) has to be replaced for this reason by 

X(l,a) = x{6 m (a)). (3.25) 

For large loops x(-) the path integral over reparametrizations in Eqs. (3.16) or (3.21) 
has a saddle point, which is denoted above by 6*(a). The saddle-point value of the 
boundary action recovers the minimal area 5' m i n [a;(-)], reproducing the exponential in 
Eq. (3.23). To say in other words, minimizing the exponential of the Douglas functional 
over reparametrizations yields the precise area-law 



K 



f ds f ds 



, [x(t(s))-x(t(s'))] 



min t(s) |e - - ^ j = e^ 5 -^. (3.26) 

We shall use this equation in Sect. 5 when calculating scattering amplitudes in QCD. 

The Gaussian fluctuations around the saddle-point 0* (a) result in an appearance of a 
pre-exponential factor, so Eq. (3.23) is modified as 

* [*(.)] lars =° ops F \VKx(-)] e -**"*[*<■)] [i + o ((KS^)- 1 )} , (3.27) 



where the pre-exponential factor F yKx{-) is contour dependent. Its calculation for a 
circle is performed in Appendix C. 

Therefore, the asymptotic area law (3.23) is recovered by Eqs. (3.12) [(3.13)] or (3.16) 
[(3.21)] modulo the pre-exponential which is not essential for large loops. 
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3.4 Remark on the area law for large loops in QCD 



In general, the Wilson loops are not observable in QCD. Observable quantities are scatter- 
ing amplitudes which are given by the sum over paths of the Wilson loops (see Eq. (2.13)). 
Nevertheless, some observables can be directly expressed through the Wilson loop of a 
certain shape. An example is the interaction potential between static quarks which is 
determined by a T x R rectangular loop for T ^> R. 

It is well-known [26, 27, 28] how this potential is calculable for the Nambu-Goto string, 
including the Luscher term. Then the most convenient choice is to parametrize the string 
world-sheet by the coordinates along the T- and i?-axes, which are conformal. The results 
described in the previous subsection allow us to answer a natural question as to how the 
linear potential can be reproduced for the unit-circle parametrization which does not obey 
in general the conformal gauge. The reparametrizing function 6* (a) is determined for the 
rectangle loop by the Schwarz-Christoffel mapping. Once again, the path integral over 
reparametrizations in Eq. (3.16) is crucial to identify the asymptote of the string disk 
amplitude with the asymptotic area-law behavior of the Wilson loop in QCD. 

4 Derivation of the Koba— Nielsen amplitudes 

Scattering amplitudes in open string theory are conventionally obtained by inserting ver- 
tex operators in the path integral over string fluctuations. Having the boundary action, 
they are represented as the path integral over x(-) with the vertex operators inserted at 
the boundary. The disk amplitude plays, therefore, the role of a generating functional for 
scattering amplitudes. 

4.1 The Fourier transformation 

It is convenient to perform a (functional) Fourier transformation 3 



is invariant under reparametrizations. This functional Fourier transformation is of the 
type as in Eq. (2.20) and transforms the position-space disk amplitude to a momentum- 
space one. 

Substituting (3.16) into Eq. (4.1) and performing the Gaussian integration, we arrive 
at the following momentum-space disk amplitude 




(4.1) 



where the exponent 




(4.2) 




m*i))-M*2))] 2 ). (4.3) 



— oo 



3 Such a transformation of the Wilson loops was first advocated in Ref. [29]. 



14 



Note that (4.3) has the same form as (3.16) only with K replaced by 1/K = 2na'. 

In front of the exponential in Eq. (4.3) there is in fact the determinant (in d — 4 
dimensions) 

-i -d/2 



detG(s(ti; 

<1*2 



= exp j-^y dtilog[G(s(ti)- 
= expj-^ J dt 1 log[G(0)]| , 



s(t 2 = h))] 



(4.4) 



which is an infinite constant. It can be regularized for the unit-circle parametrization by 
modifying the Green function (3.7): 



1 1 x ^ b m 

Gb(cr) = log (l + b 2 — 2b cos a) = — } — cosmcr , 

2"7T 71 ^-^ m 

m=l 



which yields 



detGH^)-^)) 



-d/2 



exp{-d7rlog[G fc (0)]} . 



(4.5) 



(4.6) 



What is important is that this regularized determinant does not depend on the reparametriz- 
ing function a (9) and therefore is an overall constant. 

The main advantage of the momentum- dependent amplitude (4.3) is that the momen- 
tum variable can be chosen to be a step function of t: 

M 

AO = ( 4 -7) 

where fj = 1 at the j-ih interval [tj-i,tj] and vanishes otherwise as in Eq. (3.14). Note 
that the stepwise discretization of x(t) itself is not possible since it would violate the 
continuity of the world-line of the string end. 

Since p(t) = Pj at the j-th interval, the exponential in Eq. (4.3) is in fact a function 
of M variables sf Sj < Sj+i- The only effect of the reparametrization (3.17) is then to 
change the values of s/s keeping their cyclic order: 



{si < ... < Sj-i < sj < . . . < s M } — > {h < ■ ■ ■ < tj-i < tj < . . . < t M } 
This is a discrete version of the transformation (3.17). 



(4.8) 



4.2 From disk amplitude to Koba-Nielsen amplitude 

The stepwise discretization of the momentum-space loop from the previous subsection 
naturally results in the Koba-Nielsen amplitudes. We present here a simple on-shell 
derivation which is pretty close to the standard derivations in string theory. 
Integrating by parts as in Eq. (3.22), we rewrite the exponent in (4.3) as 



+oo 



dsids 2 
(«i - s 2 ) 2 



[p(t( Sl )) - p(t(s 2 ))} 2 = - / dt 1 dt 2 #(t 1 Mt 2 )log[ S (t 1 )-s(t 2 )] 2 . (4.9) 
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For the step functions (4.7), when 
we find 

V(A Pl ,...,A PM ) = l[(l-br' A ^ I Hd Sj l[\s k - Sl \ a ' A ^. (4.11) 

* s 1 <...<Sj- 1 <s j <...<s M J fc ^' 

A few comments concerning Eq. (4.11) are in order. The factor in front of the integral 
comes from the multipliers with k — I. They are made finite by a regularization of the 
type in Eq. (4.5). The integration over s/s results from the integration over s(tj) in the 
path integral over reparametrizations in Eq. (4.3), while the integration over s(t) at the 
intermediate points is naively ignored. This integration is in the spirit of the integration 
over the positions of vertex operators in string theory. Setting a'Apj: = 1 (i.e. imposing 
the tachyonic on-shell condition) and fixing the PSX(2; R) invariance, appearing after this 
setting, in the standard way, we obtain the on-shell tachyonic amplitude in the Koba- 
Nielsen variables. 



4.3 Projective-invariant off-shell amplitudes 

As is well-known, the amplitude (4.11) is invariant under a transformation from the 
PSL(2; M) projective linear group 

s ^ s / = ^±^ ad-bc=l (4.12) 
cs + d v ' 

only when a'Apj: = 1, i.e. only for tachyonic amplitudes. On the contrary we might 
expect that an M-particle amplitude, generated by the Fourier transformation of the 
reparametrization-invariant momentum-space disk amplitude (4.3), should be projective- 
invariant because the projective group is a subgroup of reparametrization transformations. 

There are two reasons why (4.3) resulted in Eq. (4.11). The first reason is the diver- 
gence of the double integral over s± and S2 in the exponent in Eq. (4.3) for s\ = S2- This 
integral is of course convergent for smooth p(t(s)), but it is divergent for the stepwise 
p(t(s)) when si and S2 lie on adjacent sides k — I ± 1. The correct procedure is to un- 
derstand this integral according to the principal-value prescription which lead us to the 
prescription to omit the adjacent sides with k — I ± 1, as is shown in Appendix D. 

If we repeat the calculation omitting the sides with k — I ± 1, then the integrations 
over si and S2 are perfectly finite resulting in 

Sk si 

A (Pk - Pi) 

ds 2 



si - s 2 ) 2 



£ I d3l I 

= -2 • A Pl log \s k - Sl \-2j2 ^ log {Sj " 8 '-- l)(8 '" +1 " Sj) (4.13) 
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which is projective invariant. 

The second reason for the loss of the projective invariance in Eq. (4.11) is that the 
integration over reparametrizations at intermediate points of a side was not taken into 
account. A rather subtle issue how to accurately integrate over s(t) at the interme- 
diate points is explained in Appendix D, where it is shown that this integration over 
reparametrizations at the intermediate points results in the following measure 

M A 
D {M) S = TT d Sj 

II. ., _ s 



for the integration over s/s. It is invariant under the projective transformation (4.12) 
and gives 



*(Api, . . . , Ap M ) 

TT „ la'Api-Ap, TT / ( s i ~~ s i-i)( s i+i ~~ s j) 



I 



n^f^ni^.i"' aaafi n 



Sl<...<Si-l<Si<...<S M 



\ S i - S i-l\ ^ j V [Sj+i-Sj-ij 



(4.15) 



For the case of 4 scalars, Eq. (4.15) reproduces the Veneziano amplitude 

i 

A(A Pl ,Ap 2 ,Ap 3 ,Ap 4 ) = J dxx-^-^l-x)-*®- 1 , (4.16) 

o 

where a(t) = a't and 

s = -(Api + Ap 2 f = -(Ap 3 + A Pi f , 

t = -{Ap 2 + Ap 3 ) 2 = -{A Pl + Ap 4 ) 2 (4.17) 

are usual Mandelstam's variables (for Euclidean metric). Here the tachyonic condition 
a'Ap 2 = 1 has not to be imposed. 

4.4 Ambiguities of the measure 

In Ref. [14] we used, instead of (4.14), another choice of the measure 

m , v 

DM 8 = f[ds i f ^i±i_g!zij (4.18) 

f=T [ s i ~ s i-l){ s i+l ~ s i) 

which is also invariant under the projective transformation (4.12). This results in the 
amplitude 



/ 



ii ii i«. - '.i" 4 * 4 * ii ( (sj v" i -t i 1 r , T ' A "'" 



Sl<...<Si_l<Si<...<S M 

(4.19) 
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These formulas are known as the Lovelace choice [30], that reproduces some off-shell 
dual amplitudes known since late 1960's (for a review see Ref. [31], references therein and 
subsequent papers [32, 33]). For the case of 4 scalars, Eq. (4.19) reproduces the Veneziano 
amplitude (4.16) with a(t) = 1 + a't. It has the same intercept of the Regge trajectory 
as the on-shell amplitude (4.11) but now the tachyonic condition a'Apj = 1 has not to 
be imposed. 

The measures (4.14) and (4.18) can be generalized as 



D (M) S = TT ds * ( £ ! ±i^^i ) (4 2Q) 

- L .- L [Si - Si-!) V S i+ i - Si J 

where a is an arbitrary constant. The measures (4.14) and (4.18) are reproduced for 
«o — and a$ = 1, respectively. The measure (4.20) is in the spirit of the Koba-Nielsen 
amplitudes and is projective invariant for an arbitrary a®. It is the measure (4.14) which 
we have obtained in the previous subsection by integrating over reparametrizations, but 
we cannot exclude that a similar procedure may also exist for other measures. 

The amplitude associated with the measure (4.20) is again the Veneziano ampli- 
tude (4.16) with a(t) = a + a't. Therefore, only an intercept of the Regge trajectory is 
sensitive to the choice of the measure provided it is projective invariant. But the change 
of the measure (4.20) becomes not essential for I /a' <C — t < s, when the integral over x 
is dominated by a saddle point. 

The Regge trajectory with a(0) ~ 1 is usually associated with the vacuum trajectory 
= Pomeron which appears in QCD from cylinder diagrams ~ 1/N 2 (with two quark 
loops). We work with planar diagrams, single quark loop and, correspondingly, with a 
quark- ant iquark Regge trajectory which has a(0) ~ 1/2 from experiment. 



tt 



5 QCD amplitudes and the area behavior of W 

In the following we ask what is the contribution from the area behavior of the Wilson 
loop to the corresponding QCD amplitudes. It is of course clear that depending on the 
kinematical situation there will be other important contributions, e.g. at large transverse 
momenta where a perturbative behavior of W is relevant. 

As discussed in the introduction it is well-known from lattice calculations that the area 
behavior, and more generally large distance results, have been found to follow the predic- 
tions from the Nambu-Goto action rather precisely even at surprisingly low distances. It 
therefore makes sense to ask what happens to the QCD amplitudes if we approximate the 
Wilson loop by the area behavior. This question will be studied in the present section. 

We now want to insert the Wilson loop given by the area behavior, taking the latter 
from the celebrated Douglas construction [25] discussed in details in Sect. 3 and Ap- 
pendix B. We write 

W = SV a j Va{r) exp (J^- &d d9'z(0)z(0') ln(l - cos[<t(0) - <t(0')])) ■ (5-1) 
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Here SV a means that the full integral should be replaced by the Saddle Point with respect 
to the reparametrizations a ((f)), where a satisfies & > and er(0) =0 and a(2ir) = 2ir. 
The integral over er(r) would then produce the exponential of the minimal area times a 
prefactor as displayed in Eq. (3.27). However, since we want to consider QCD, then only 
the minimal area factor is pertinent, as we know from lattice gauge calculations. In QCD 
there are sub leading factors (perimeter term, Liischer term, ...), but they are probably 
not given by the prefactor coming from the saddle point in (5.1). Therefore the SV a 
is restricted to include only the leading saddle-point contribution with no multiplicative 
prefactors. 

In the following we want to insert Eq. (5.1) in the basic formula (2.20) with the 
goal of doing the z— integral. If we perform the Douglas saddle point SV a first, then 
the resulting minimizing functions a will depend on the curves C , and hence also on 
parameters entering in z. As an example, in Appendix B we have discussed the elliptic 
case, where the minimizing a depends on the ratio of the lengths of the two axes, and these 
lengths also enter in the expression for z. Hence, following this procedure the z— integral 
cannot be performed as a simple Gaussian integral. Therefore, in applying Eq. (2.20) we 
shall assume that the SV a operation commutes with the z— integral entering in (2.20), 

J Vk{r) J Vz{t) SV a J Va(r) = J Vk{r) SV a J Va(r) J Vz(r) , (5.2) 

so that finding the saddle point can wait until the k— integral is to be performed. 4 The 
validity of Eq. (5.2) is plausible a posteriori, since the a— dependence in all situations 
encountered turns out to occur entirely in the logarithm ln(l — cos(er(r) — cr(r'))), as we 
shall see in the following. 

With this assumption we can now perform the Gaussian z— integration in the basic 
formula (2.20) to obtain 

roc pT M ~ 2 pr m p p 

G(A Pl ,...,A PM ) = j dTe- mT j dTM-iHj d n J Vk{r) SV a J Va(r) 

X Sp P g ^/o Td ^/o Td ^^M+PW)-(M T OW))ln(l-cos( ( T(r)- ( T(r')))-i/ ( rdr 7 (r).fc(r) > (53) 

Here it should be emphasized that the dots in Eq. (5.1) can be moved to the logarithm 
by partial integrations, and hence the integration over z can be formulated such that it 
involves only z itself and not i. We mention that a similar simplification does not occur 
in Eq. (2.23) for scalar quarks, due to the occurrence of z 2 . 

Inserting the stepwise p(r), regularizing the integral in the exponent by the principal- 



4 It worth noticing that the SV„ operation is in fact nothing but taking the classical limit h — > 0. 
If the dependence on Planck's constant h is restored, the exponents in both the coordinate-space disk 
amplitude (3.16) and the momentum-space disk amplitude (4.3) are divided by h because it enters the 
exponent of the Fourier transformation. 
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value prescription and using the formula of the type of Eq. (4.13), this becomes 



M-l M 

G(A Pl ,...,Ap M )oc J] / dfa n 

1 { 3=1 

( 



x exp 



1 M 

— ApiA Pj ln(l - cos(0i - fa)) 



sin[(0 J+1 - fa)/ 2] sm[{fa - <ftj-i)/2] 
sin[(0 i+ i - 0j_i)/2] sin 2 (fa /2) 

\ 



Ap2/47rii- 



AtxK 



\ 



*>j = l 
»*7 



/ 



00 = 0, 0M = 27T, 

where the "kernel" K, is given by 



/C(0i, . . . , 0m-i; Api, . . . , Ap M ), 

(5.4) 



/C 



Pife(0) SV a J Va((j)) 



x exp 



x exp 



AivK 



2?r 



2tt 



d6 / d& k(9)k(9') ln(l - cos(<t(0) - <r(0'))) 



d0A;(0) ln(l - cos(0i - <r(0)) 



/ 



x / dr r 



M— 1 „ -mr 



sp P exp 



2tt 7o 



d</> 7 (</W) 



(5.5) 



Minimizing the kernel /C with respect to a gives the requirement 



2tt 



| d0' fc(0)fc(0') cot ^( g )-^( y ) + 2 £ A Pl fc(0) cot = 0. (5.6) 

o * 

For a given k this determines <7*(0) or, alternatively, the minimizing reparametrization 
0*(er). Using er* it follows by our assumption (5.2) that only the minimal areas are included 
as contributions to the amplitude. 

The main result (5.4) reveals the interesting appearance of poles of the Veneziano 
type. These occur when two fas coincide, fa — > <pj. The kernel K, does not in general 
have any pole singularities dependent on the momenta when fa — > fa, so the momentum 
dependent Veneziano-type poles cannot be canceled by any contribution from K. 

The factor 



exp 



— j^: A Pi A Pj H 1 ~ cos (0i - <t>j)) J 



(5.7) 



is independent of the reparametrizations involved in the Douglas construction. This is 
due to the fact that p is a simple sum of delta functions. The above factor is therefore 
universal. 
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After a transformation of the variables fa, Eq. (5.7) becomes similar to the Koba- 
Nielsen representation of the M-point function. To see this use 5 

In [2(1 - cos(0i - fa))} = 2m|2sin((0 i - fa)/2)\ = 2 In | Si - Sj\ + . . . , s t = -cot(&/2), 

(5.8) 

where we left out terms that vanish in the sums occurring in Eq. (5.7) due to energy- 
momentum conservation. The Sj's occur like the variables in the Koba-Nielsen formula. 6 
Thus, in QCD the angular variables, related to the proper-time variables by fa = I^Tijr, 
play the role of the Koba-Nielsen variables. 

We thus see that the usual dual model poles are present in the QCD amplitude (5.4). 
This is indeed to be expected on intuitive grounds, because the area behavior of the 
Wilson loop can be interpreted in a string framework as arising from the rotating stick 
with a well-known Regge type spectrum. 7 It is therefore quite satisfactory that this result 
also occurs in our general QCD formula when the area law is imposed. 

It must be emphasized that in QCD the area law is valid only for large areas, i.e. for 
large distances. Translating this to momentum space, we need large momenta. Therefore, 
the low lying dual model spectrum is not relevant in QCD. Thus, the tachyon is also of 
no relevance, as should be the case. 

Going back to Eq. (5.4), it should be noted that this formula does not correspond to a 
simple dual amplitude, because of the additional factor K. Thus, although the spectrum 
is quite stringy, the amplitude is more complicated than in the standard dual models. 
Presumably, this is not too surprising. 

As is mentioned above, the tachyon does not occur due to the stringy correspondence 
large distances=large momenta (in contrast to perturbation theory). Also, when the 
area behavior leads to extremely small contributions, as is the case for large transverse 
momenta in string amplitudes (for the 4-point function — t ~ s) , these are not important 
relative to perturbative contributions 8 , which would then dominate the Wilson loop. This 
is how the exponential falloff of the 4-particle amplitude with large —t ~ s, which is 
unavoidable in string theory [36], does not happen in our consideration. 

Now let us ask if it is possible (at large momenta) to obtain something like a standard 
dual model. This depends on the extra factor K, in Eq. (5.4). The last factor in the 
integral over k in (5.5) suggests that we substitute 

k -> k/r. (5.9) 

Then we see that the first two factors in the definition of JC become close to one if r is large. 
For the number of external mesons M large, the r— integral in Eq. (5.5) is dominated by 
large 

T = (M-l)/m. (5.10) 

5 Here and below the value of 4>m can be chosen arbitrary, respecting the cyclic symmetry. 
6 The differentials transform as d<f)i = 2dsj/(l + sf). 

7 Another derivation of the Regge spectrum from the are law was given in Ref. [34] . 
8 We would like to emphasize once again that we are dealing in the large N limit with the quark- 
antiquark Regge trajectory, whose perturbative QCD calculation was pioneered in Ref. [35]. 
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It is seen from this formula that the values of r, dominating the r— integral in Eq. (5.5), 
are also large for small values of m, but we consider such a limit to be rather formal 
because m has the meaning of a constituent quark mass in QCD, which is about hundred 
MeV from experiment even for very small bare masses of up and down quarks because of 
the spontaneous breaking of chiral symmetry. 

When the large values of r (5.10) dominate, the integrand of the integral over a (fa) in 
Eq. (5.5) does not depend on er and /C degenerates into 



f 7W^ TT sin(<fc + i/2) sin(<fc/2) tt 1 ( , in 

a(<t>i)=(t>i 

modulo a constant which does not depend on the fa's. The appearance of this factor is 
due, shortly speaking, to the specifics of the Douglas minimization for stepwise functions 
p(0(<f>)) in contrast to smooth functions. The points fa's, where the function has dis- 
continuities, are irregular points from the point of view of the minimization because the 
minimizing function 9* (fa) has to satisfy 9*(fa) = fa. We can still perform a reparametriza- 
tion <f> — > 9 (fa) at intermediate points G (fa, fa+i). For such a reparametrization p(9(fa)) 
moves along the step but its value remains unchanged. Therefore, this is a zero mode 
in Douglas' minimization and we have to integrate over these zero modes. For smooth 
functions, 6*(fa) was just fixed. The integration over the zero modes is exactly the same as 
the integration over reparametrizations at intermediate points described in Subsect. 4.3 
and Appendix D. The result of this integration at the interval (fa, fa+i) is given by the 
i-th multiplier in Eq. (5.11), while the product runs over the labels of intervals. 
The scattering amplitude then takes the form 

M-l S< + J , M , n 11 in Ap?/27rK 



G(A Pl ,...,Ap M ) oc I]/ 

1=1 „ 



dsi 



1 1=1 



1 




s i+l ~~ s i\ 




\ s i+l ~ s i 






Si+l - 


' Si-i 





1 M 

X 6XP iz^K ^ A P iA Pi ln \ s i- s i\) ( 5 - 12 ) 



<>3 = 1 



which looks similar to the Koba-Nielsen amplitude. It should be kept in mind that this 
expression is only valid if the area-law behavior of the Wilson loop dominates over other 
contributions to W , and if the number of external particles M is large and/or m is small 
compared to \fK as is already mentioned. 

In Appendix E we have discussed this expression in much more details. In particular 
we have shown how the Veneziano amplitude with the usual Regge behavior follows. 



6 Large number of external particles WL/SA duality 

The reader may have noted that the area behaved Wilson loop and the large M amplitude 
before the integrations over the fa's are done look very similar (see Eqs. (5.1) and (5.12)), 
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except that the variables z and p are somehow interchanged. This reminds us about 
the Wilson-loop/scattering-amplitude duality mentioned in the introduction, which was 
found in SYM perturbation theory, see Refs. [8] and [9]. 
This similarity can be made explicit by the substitution 

i 

where the O's are the Heaviside step functions. From this expression we have 

m = ^ E a p* ^ - = it ? Apt W) _ 9i) ' (6 - 2) 

i i 

where 9 = 9((f>) is some reparametrization with 9(<f>i) = 9i. If this equation is inserted in 
the area behaved Wilson loop (5.1), it reproduces the multiparticle amplitude in (5.12) 
when integrated over the 0j. The spatial variable on the left-hand side of Eq. (6.1) is on the 
right-hand side composed of a number of constant vectors given by the momenta, which is 
rather analogous to the duality (1.2) mentioned in the introduction for the M = 4 SYM. 
A difference is that here we need to integrate over all trajectories considered as functions 

of (pi. 

Equation (6.1) reemphasizes the often mentioned fact that for strings the large z limit 
is equivalent to the large p limit. The surprising interchange of space and momentum is 
clearly a stringy effect. 

It is easy to check that the expression (6.1) actually satisfies Douglas' variational 
condition (B.5) or, alternatively, Eq. (B.14). To see this, let us make a variation 

(f)(9) = &(0) + <ty(0), (6.3) 

with 5(f) small and 0* giving the minimal Douglas functional, as is explained in Ap- 
pendix B. We have the conditions 

z(9,((f) i )) = Xi , (6.4) 

since the curve passes through the points Xj. Therefore in the variation we need the 
conditions 

= (6.5) 
for all %. The first variational derivative of the Douglas functional is proportional to 

J d9 J d9' 5(f)(9)mH0') cot M^-M?) ^ (6 6) 



Inserting Eq. (6.2), the expression (6.6) becomes 

E A p* A Pi cot ( 6 - 7 ) 

ij 
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where we get contributions only from the fixed points X{. Because of Eq. (6.5) we see that 
the first variational derivative (6.7) vanishes. Hence the duality relation (6.1) actually 
satisfies Douglas' variational principle and in this sense the curves (6.1) represent the 
dominant trajectories (the "master trajectories") in phase space when the momenta are 
given. The special trajectory (6.1) is a kind of zero mode solution of the Douglas varia- 
tional problem, since it does not actually determine the function 0*, which is irrelevant in 
the case of stepwise constant momenta. We have already discussed this specifics of step 
functions in Sect. 5. 

To sum up, the result is 

M-l &+ 1 m . 2( , /t) \ 

G(A Pl ,...,A PM ) oc J] / d ^II • rrT Iwoi 

x w(z(<t>) - ^ $>9(0 ~ ^)e>(^+i - 0)) , (6-8) 

i 

where W(z(<j>)) is the Wilson loop as a function of the boundary curve z(<f>). Except for 
the integrations this is similar to the supersymmetric case (1.2) discussed in [8] and [9]. 

This similarity can be made even more explicit by noting that Eq. (6.1) means that in 
the interval fa < <fi < fa + i the vector z{fa) is equal to the vector Pi/K. Remembering that 
z equals X{ and Xj+i for the parameter values fa and fa+i, respectively. This amounts to 

Api = K - Xi), for fa < <j) < fa +1 , (6.9) 

in conformity with Eq. (1.2). 

We again emphasize that all the above results are valid only when the number of 
external particles M is large. In the general case of a smaller number of produced particles 
the situation is much more complex. Presumably there will be important fluctuations 
around the master trajectory giving contributions to the kernel K. 

7 Conclusions 

We have found a relation between the meson scattering amplitudes and the Wilson loop 
for large iV QCD. We then investigated the behavior of an area behaved W when a func- 
tional Fourier transform was performed, leading to the Veneziano multiparticle amplitude 
in the Koba-Nielsen formulation. This turned out to be very useful when we inserted an 
area approximation for the Wilson loop in our general expression for the large N QCD 
amplitude. The result is a convolution integral, with the well-known Koba-Nielsen inte- 
grand convoluted with a kernel /C. The usual poles always occur because they cannot be 
prevented by the kernel. 

For a very large number of external particles JC becomes essentially a constant, and 
hence the Veneziano multiparticle amplitude appears. It then turns out that there exists 
a nice duality between Wilson loops and scattering amplitudes, somewhat similar to the 
supersymmetric case. 
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Although these phenomena are valid for large N, one might hope that something 
similar occurs for N = 3, in which case this would be observable at LHC, where a huge 
number of particles are produced. So we hope to see at least some tracks of the Veneziano 
amplitude in the collider data! 
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Appendix A An example of the use of the "variable" 

gamma matrix 

Here we shall give an example of how the gamma matrix 7 At (t) operates in the case where 
we consider the quantity 



F = sp P exp 



/ di p„(t)7„(*) = § p n ex p - - ( a - x ) 



where we took the momenta p to be stepwise constant. We can now use 

fa- 

e VjMj = cosMAtj + sin MAtj, Atj =t j -t j - 1 . (A. 2) 

Here we took p 2 = — M 2 for all j, thus assuming that all the external mesons have the 
same mass. Then 

f = sp n ( cos ma ^ + 1 m sin MAtj ) ■ ^ A ' 3 ^ 

j=n ^ ' 

For the two point case n = 2 we get 



Also 



F n=2 = 4 ( cos MAti cos MAt 2 - sin MAt x sin MAt 2 ) . (A.4) 

17 A( P^ 2 P^ PlP* \ / A r\ 

F n=3 = 4[c lC2 c 3 - jpS lS2Cl - j^s 1 c 2 s 3 - ^c lS2 s 3 ) , (A.5) 



and 



j? J P3P4 P2P4 P3P2 PlPA 

F n=4 = 4^CiC 2 C 3 C 4 - ^^CiC 2 S 3 S 4 - ^-CiS 2 C 3 S 4 - ^-CiS 2 S 3 C4 - ^-sic 2 c 3 s 4 

!> ' P3 SlC 2 S 3 C 4 - ^^SlS 2 C 3 C 4 + y^((plP 2 )(P3P4) - (PlP3)(P2P4) 



M 2 M 2 1 ' ° * M 2 

+ (pm)(p 2 P3))sis 2 s 3 s 4 y (A. 6) 

Here q = cos MAtj and Sj = sin MAtj. 

If the momenta are not stepwise constant, we can use the above procedure with n — > oo, 
if we divide the interval from to r into intervals r/n and take the limit at the end. 
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Appendix B Douglas' approach to the minimal area 



The Douglas algorithm [25] for finding the area of the minimal surface bounded by a 
closed contour C which is parametrized by the function x^a) is based on minimizing the 
boundary functional 9 



2tt 2tt 



A lxm = -1 [do / d yMWh!Mr (B.l) 

L WJ 8?r J J 1 - cos(a -a') K ' 



with respect to the reparametrizing functions 9(a) (c\9(a)/da > 0). The numeric value of 
A for the given C depends on the choice of 0(a) and in general 

A[x(9)] > S min (C) (B.2) 

while the equality is reached for certain function 9(a) = 9^(a) which provides the minimum 
of A. The function 9*(a) is of course contour-dependent. 
To prove the fact that 

A[x(9*)] = S min (C) , (B.3) 

we reconstruct -X" M (r, a), describing the surface, in the interior of the unit disk r < 1 
from the boundary value X^(l, a) = x^ (9*(a)) by the Poisson formula and note that thus 
constructed X M automatically obeys conformal gauge 

f¥L — - o r 2 — — - — — (B 4) 

dr da dr dr da da 

Therefore, the Nambu-Goto action coincides for this configuration with the quadratic 
action and the boundary action in Polyakov string theory coincides with the area. 

By varying A[9] with respect to 9(a) at given C, we get the following equation for 
0*(<7): 

[ djl x(9*(a))-[x(9*(a))-x(9 if (ai))] _ Q 

J 1 - cos(a - a') " 1 ' ' 

o 

which can be written in several equivalent forms. 

To illustrate how Eq. (B.5) can be used to determine 9*(a), let us consider the case of 
a plane contour, when the problem can be solved by a conformal map, and concentrate 
on the case of an ellipse 

X\ = a cos 9(a) , x 2 = b sin 9(a) . (B.6) 
Then Eq. (B.5) takes the form 

J da sin [9,(a + a) - 9,(a)] _ ^ J sin %{a + a) + 9*(a)] - sin [20»(<r)] 
J 1 — cos a J 1 — cos a 



3 For a modern review see also Ref. [37], Appendix H. 
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with 



a 2 -b 2 



a 2 + b 2 ' 

For the simplest case of a circle (a = 6), the coordinates r, a are conformal so that 



(B-8) 



#*(<t) = a | circle 



(B.9) 



The right-hand side of Eq. (B.7) vanishes for e = 0, while the left-hand side also vanishes 
for 0*(cr) given by Eq. (B.9). 

For e 7^ the following ansatz passes through Eq. (B.7): 



9*(cr) = a + ^ c n sm(2na) ellipse 



(B.10) 



n>l 



after which a set of algebraic equations relating c n 's emerges. For small deforming pa- 
rameter e, its iterative solution to order C(e 7 ) found by Mathematica is 



Cl 

c 2 
c 3 
c 4 
c 5 
c 6 



e - e 3 /4 + e 5 /8 + 0(e 7 ) 
3e 2 /4 - 5e 4 /8 + 25e 6 /64 + £>(e 7 ) 
5e 3 /6 - 5e 5 /4 + 0(e 7 ) 
35e 4 /32 - 77e 6 /32 + £>(e 7 ) 
63e 5 /40 + C(e 7 ) 
77e 6 /32 + C(e 7 ) . 



The minimal area 



*Smin = Kab = 7TG[/ 



1-6 



(B.H) 
(B.12) 



is of course reproduced to this order by substituting (B.10), (B.ll) into Eq. (B.l). 
Integrating by parts, we can also rewrite the Douglas functional (B.l) as 



2n 



2tt 



A = -i- J d0i J d9 2 ±(0!) ■ x{9 2 ) In (1 - cos [<r(0i) - a{9 2 )\) . (B.13) 



Its variation with respect to er(0) results in the equation 



2tt 



/ 



da £(#) • x(a) cot 



cr*(0) - a* (a) 



(B.14) 



which is nothing but the Douglas original equation from the Abstract of 1927. 

By solving this equation, we obtain the function a*(9) which is inverse to 9*(cr), given 
by Eq. (B.10), and reads 



&*(9) = 9 + d n sin(2n$) ellipse 



:b.i5) 



n>l 
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with 



di = -e + 0(e 7 ) 

d 2 = e 2 /4 + e 4 /8 + e 6 /16 + 0(e 7 ) 

d 3 = -e 3 /12-e 5 /16 + C(e 7 ) 

d 4 = e 4 /32 + e 6 /32 + C(e 7 ) 

d 5 = _ e 5/80 + 0(e 7 ) 

d e = e 6 /192 + C(e 7 ). (B.16) 
To proceed further, it is convenient to use a standard ellipse of the area it with 



Then 

xi(a) = a cos 9(a) = /j n cos(2n — l)a 

n>l 

x 2 (a) = bsin6(a) = v n sin(2n - l)tr (B.18) 

n>l 

and from Eq. (B.ll) fi n = v n to order 0(e 7 ). 

This can be understood if we continue the boundary coordinates (B.18) inside the unit 
circle as 

Xi(r,<r) = ^/i n r 2n ~ 1 cos(2n- l)a 

n>l 

X 2 (r,a) = ^z/ n r 2n - 1 sin(2n- l)a . (B.19) 



n>l 



It can be then explicitly verified that these coordinates obeys the conformal gauge (B.4) 
for fi n = v n . 

Introducing the analytic function 

M(z)=J2^z 2n - 1 , (B.20) 

n>l 

that describes a conformal map of a unit disk onto the interior of the standard ellipse, we 
get the following algebraic equation 

acos6(a) +ib sin 6(a) =M(e ia ) (B.21) 



or 



e = M(e*) + y /M*(e*)-(a*-V) 

a + b 1 ' ' 
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which determines 9(a) for the given Ai(z). For the ellipse, 9 plays the role of an angular 
variable in the parametrization (B.6), while the function a*(9) relates it to the variable a 
inherited from the conformal coordinates r, a obeying Eq. (B.4). 

An analytic function that describes the conformal map of a unit disk onto the interior 
of an ellipse was found by Schwarz in 1869 [38] and gives 



where 



and 



M(z) = y/a 2 -b 2 sin 



F(z;s) 



7T 



F 



2K(s) Vv^ 



ax 



y/(l-X 2 )(l- S 2 X 2 ) 



1 

K(s)=F(l;s) = [ dX 
J J(l-x 2 )(l 



^(l-x 2 )(l-s 2 x 2 ) 



(B.23) 



(B.24) 



(B.25) 



are, respectively, the incomplete and complete elliptic integrals of the first kind. 10 The 
parameter s is related to e by 



a + b nK{VT^) 
log = 2^(s) = 



a — b 



2K(s) 



(B.26) 



so that iteratively 



8 = 2e - -e 3 + -e 5 - -e 7 + —e 9 + 0(e n ) . 
2 2 32 128 v ; 

Substituting the function (B.23) into Eq. (B.21), we obtain the equation 

e iCT 



(B.27) 



(B.28) 



whose imaginary part reproduces Eq. (B.26) in view of the important identity 



; s 



\k (yi - s 2 ) . 



(B.29) 



The real part of Eq. (B.28) determines the function 9(a). 

In order to obtain it, it is convenient first to rewrite Eq. (B.23) as 



M(z) = Va 2 — 6 2 cosh 



arccosh(z/ v / s) 



7T 



2K(s) 



dX 



o 



y/l — s 2 cosh 2 A 



(B.30) 



10 We use the notations for elliptic integrals from Wikipcdia. They are related to those of Mathematica 
as F (z; s) = EllipticF[ArcSin[z], s 2 ] and K{s) = EllipticK[s 2 ]. 
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This form is more suitable for the case when the first argument of F in Eq. (B.23) is large 
as it is for small e and correspondingly for small s, related by Eq. (B.27). Substituting 
the function (B.30) into Eq. (B.21), we obtain the equation 

arccosh( c 1<T / \/~s) 

* + * = ^ f , ^ (B.31) 

o v 

which is equivalent to Eq. (B.28). 

Differentiating Eq. (B.31), we obtain 

tft \ = 71 1 = 71 1 fB32) 

1 ' 2K(s) VI + s 2 - 2s cos 2a 2K(s) ^(1 - s) 2 + 4s sin 2 a 1 ' ' 

which yields 

„/ n t f da' 7T 1 „ / 2iv/i \ .„ , 

0(a) = —— — -_ = -F sin a ; -. v x B.33 

K> 2K(s)J ^(1-5)2 + 4s sinV 2K(s) (1 - s) \ ' (I - s) ) 1 ; 
o 

which is again an elliptic integral of the first kind. This reproduces iteratively Eq. (B.ll) 
provided Eq. (B.27) is satisfied. 

The solution (B.33) obeys the properties required for the Douglas minimizing function 
9*(a). It is seen from Eq. (B.32) that 01(a) is positive and finite for s < 1, as is required 
for a reparametrization. For s — > we have 0*(<r) — > 1 as it should for a circle. In the 
limit s — > 1 when K(l) = oo we have b — > 0, so the ellipse collapses. Then #*(cr) vanishes 
everywhere except for = 0, n where it becomes infinite: 

0i (ff) ^° 7T(J((T) + 7T<5(<7 - tt) , (B.34) 

while 



0»(<r) U 7r9(ff) + 7T0( ( j-7r) (B.35) 



is stepwise. 



Appendix C The pre-exponential in Eq. (3.27) 

To calculate the pre-exponential in Eq. (3.27), we substitute 

9(a) = 9*(a) + f3(a) (C.l) 

or 

a(9) =a*(9)+P(9), (C.2) 

where 

(3(0) = (3(2n) = (C.3) 
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and expand the exponent in Eq. (3.12) (or Eq. (3.13)) to the quadratic order in (3. The 
linear in (3 term vanishes because 6* (a) is the extremum while the quadratic part reads 

S*W)\ = y / dM0 2 f(0i)f(0 2 ) G- 1 - <r,(0 2 )) [/3(0i)/3(0 2 ) - /? 2 (#i)] . (C.4) 

The function (3(9) has to obey 

P'(6) > -a'M (C.5) 

for the derivative of the reparametrizing function to be positive. This is always satisfied if 
(3 is small and smooth enough. Therefore, the measure for the Gaussian path integration 
over 13(6) is the usual one for smooth functions (3(6) but we shall see subtleties for the 
functions with large derivative. 

In order to calculate the pre-exponential in Eq. (3.27), we need to do the Gaussian 
integral 

J 2 = Jv(3(6)e~ S2[m] (C.6) 

with S 2 [(3(6)\ given by Eq. (C.4). 
For a circle of the radius R, when 

x 1 (9) = Rcos9, x 2 (6) = Rsin6, x 3 (6) = x 4 (6) = (C.7) 

and cr*(9) = 6 according to Eq. (B.9), we have 

KR2 ^ ,„ ,„ COfi(6 1 ~6 2 



\ " d6M 1 C ^T6 ^6 \ [MM 6 *) - ■ (C.l 



Ait 

It is now possible to calculate the path integral over (3(6) by the mode expansion 

oo 

(3(6) = a + ^ (a n cos nd + b n sin nd) , (G9) 

n=l 

where a = — °« to obey the boundary condition (C.3) and a n , b n have to satisfy 

oo 

^^n(—a n sin n6 + b n cos nd) > — 1 (C.10) 



for the restriction (C.5) to be fulfilled. Inserting Eq. (C.9) into Eq. (C.8), we find 

ttKR 2 



2 



E^-^K + ^n)- (Cll) 



n=l 



A consequence of Eq. (C.ll) is that a\ and b\ are zero modes. They are, however, 
restricted by (C.10) as —1 < a 1: bi < 1, so we believe that the integrals over a x and &i 
simply give a constant. The integrals over nonzero modes are Gaussian and if we were 
not take into account the restriction (C.10), the result would be 

oo „ oo 

h oc J] / da ^n e 52 = n [KR 2 (n - 1)/2] _1 . (C.12) 

n=2 n=2 
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The infinite product in Eq. (C.12) can be calculated in the standard way, using a 
regularization via the ^-function: 

oo 

n=l 

oo 

Y[n = \Fhx. (C.13) 

n=l 

We then obtain from Eq. (C.12) 

h oc (KR 2 ) 1 ' 2 (C.14) 

which fixes the KR 2 dependence of the pre-exponential in Eq. (3.27). 

Since the typical values of f3, which are essential in the path integral over (3 in Eq. (C.6), 
are (3 ~ l/y/~KR, i.e. small for \f~KR » 1, the higher terms of an expansion of A[9*(a) + 
0\ in (3 are suppressed at large \f~KR. The loop expansion goes in the parameter 1 / KR 2 
and only one loop contributes with the given accuracy. We thus reproduce the behavior 
of the type in Eq. (3.27). 

We can now ask the question whether or not the typical values of a n 's and 6 n 's, which 
are essential in the integral over the nonzero modes in Eq. (C.12), obey the restriction 
(CIO). We estimate these values as 

a n ~b n ~ (nKR 2 )- 1/2 (C.15) 

which are small in accordance with the standard wisdom that high modes are not essential 
in the Gaussian path integral. However, the restriction (C.10) can be satisfied by a single 
mode only if 

n < KR 2 . (C.16) 

We assume that this type of the restriction on the number of modes becomes unimportant 
for KR 2 -> oo. 

Appendix D Path integrals over Diff(M) 

The measure for integrating over reparametrizations is determined by the metric [21] 

+oo 

\\6t\\ 2 = J dst'(s) [5t(s)} 2 (D.l) 

— oo 

which is invariant under reparametrizations. It differs from the usual one 

+oo 

INI 2 = Jds [5t(s)] 2 (D.2) 

— oo 

by the presence of t'(s) = dt(s)/ds in the integrand. 
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An explicit representation of the measure £>diff t(s) for the integration over reparametriza- 
tions is given by Eq. (3.18). As distinct from the usual measure Vt{s): 



+00 L 



i Vt(s) ■ ■ ■ = lim I J[ dtj ■ ■ ■ , (D.3) 

J L^oo J 

— oo J 

the integrals in Eq. (3.18) are ordered. 

The measure (3.15) is the invariant measure on the group Dif f{S 1 ) of reparametriza- 
tions (diffeomorphisms) of a circle. Analogously the measure (3.18) is the invariant mea- 
sure on the group Diff(M) of reparametrizations (diffeomorphisms) of the real axis. They 
are defined in the way to be invariant under the PSL(2; M) projective transformation at 
very small but finite discretization spacings £j = Sj — Sj_i ~ 1/L. This guarantees the 
invariance of the measure under reparametrizations in the limit L — > 00. 

If the integrand of the path integral over reparametrizations is not a functional but 
a function of, say, only t(sj) = tj and t(sk) = t k with j > k, we can integrate over 
intermediate values Vs with k < i < j. Analogously, for a function of M variables we 
integrate over tj's inside M intervals. The result of such an integration is a function of 
the remaining M variables t^, . . . , t iM , which should be covariant under the projective 
transformation. 

The simplest integral we meet is of the type 

U+i 

I (U+i — U){ti — U-i) ^ ^ 

ti-1 

which is logarithmically divergent at the upper and lower limits of the integration. We 
regularize it by introducing a small 5 as 

ti+i 



/& 2 
dtj -, — ; — 77 —. — t = . (D-5) 

(t i+1 -uy-'iu-ti-!) 1 - 6 



lim 

5 



We can now perform a heuristic derivation of the measure (4.14), which is invariant 
under the projective transformations, using the following formula 

J (tK-uy-^iu-toy-^ (tK-to) 1 -^ v ; 

to 

Equation (D.6) is an analogue of the well-known formula 

+00 



e -(t K -t i ) 2 /2u 1 e -(ti-t ) 2 /2u 2 e -(t K -t Q y/2(iy 1 +u 2 ) 



— LXJ 

which is used for calculations of path integrals with the usual Wiener measure. 



(D.7) 
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Choosing in Eq. (D.6) v\ = S\5 and u 2 = 628 and repeatedly integrating over interme- 
diate points ti, . . . , tx-i, we obtain 



K-l u+1 
limJJ / dt, 



i=i 



£1 



2^i=j £ j 



,s K - s 



to 



(u+i - uy-^ (*! - t y-^ (t K - 1 ) (t K - 1 ) ' 



(D.8) 



Equation (4.14) can be now derived fixing certain M values of tj's and repeatedly inte- 
grating over the intermediate points. 

In order to justify Eq. (4.15), we again introduce L infinitesimal intervals E{ and 
consider a step function p(t(s)) which has discontinuities only at finite number M of 
points t K . (J — 1, ... , M), so Api = at all other points which we call the intermediate 
points. To emulate the principal- value integral in the exponent in Eq. (4.3), we omit the 
terms when Si and s 2 he at two adjacent infinitesimal intervals. Similarly to Eq. (4.13) 
we get 



tfc+i t(+i 



d«i / ds 



k,j=l , 
k^l±l tk 



= lim 

L— >oo 



(Pk ~ Pi) 
(si - S 2 ) 



-2 A Pj ■ A Pi lo § \tj~ti\-2^2 Apj log 



3,1=1 



= -2 



M 



M 



a Pk 3 ■ Ap Kl log \t Kj -t Kl \ + Y^ A Pk 3 1o S 



(t j+ i - - 
(tj+i - tj-i) 

(tK j+1 -t K] _^){t K] -t Kj _ lt 



3,1=1 



K 



3+1 



tKi-i, 



(D.9) 



where we have substituted t,-'s for s/s and taken into account that Api = for the 
intermediate points with i 7^ Kj. 

We still have to insert (D.9) in the exponential and to integrate over the intermediate 
points. Let us consider a piece from the point K\ to the point K 3 given by the ordered 
integral of the type 



lim 

8^0 



K 2 -l 



\tKi — 



n ^ 



X 



i=K 1 +l 
e K 2 



\t ■ — t ■ 1 1 l 1 -^ 5 



1^2-1 _ ^K 2 \\^K 2 _ tK 2 +l I 



I** 



2-1 



tK 



2 + 1 I 



K 3 -l 



\tK 2 -t K2+l \ 1 ^* 



n 



5,; <5 



j=x 2 +i 



/• — /•II l 1 -^ 5 



(D.10) 



The integrals over all intermediate points, except for the ones with % = K 2 — 1 and 
% = K 2 + 1, are the same as in Eq. (D.8) and are easily doable. The two remaining 
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integrals are 



;d.io) 



lim 

<5^0 



dt 



Si 



£k 2 -i$ 



K 2 +r 



t Kl - t K2 -i\ l - £lS |*tf 2 -i - JtfJ 1 - 4 **- 1 * 



x 



: — tK 2 \\tK 2 ~ L K 2 +: 

1**2-1 ~~ *X 2 +1 | 



*; 



where 



if 2-2 

£\ — ^ £ » ) 

j=*i 



- «'^P 2 K 2 



e K2 8 



\t K2 - t^+il 1 - 6 ^ 6 IWi - ^ 3 l 1_£l5 ' 

(D.ll) 



ifs-1 

^2 = £i 
i=K 2 +l 



(D.12) 



For nonvanishing Ap|- 2 the integrals in Eq. (D.ll) differ from the one in Eq. (D.8). 
However, they are also easily calculable for 5^0 when only the domains (** 3 — ** 2 +i) ~* 
and (** 2 -i — **J — > contribute. We get finally 



(D.ll) oc 



1 


**1 — ** 2 


**2 ~~ ^3 


**1 — ** 2 


^1 " 


- ** 3 



** 2 — *Ks 



(D.13) 



thus proving Eq. (4.15). 



Appendix E Regge behavior of QCD scattering am- 
plitude (5.12) 

We show in this appendix that the scattering amplitude (5.12) has asymptotic Regge 
behavior and reduces to the Veneziano amplitude for M — 4. 

The integrand in Eq. (5.12) is the same as in Eq. (4.15) except for the additional 
factors 1/(1 + sf) in the measure. For the 4-point function we obtain 



G 4 



S4 S3 S 2 S4 S3 +00 

ds 3 J ds 2 J dsi + J ds 3 J ds 2 J ds 1 

-oo — oo — oo — oo — OO S4 

oo s 2 +oo s 3 s 2 



+ 



J ds 3 J ds 2 J dsi + J ds 3 J ds 2 J ds : 



— OO S4 S4 

1 



S4 S4 S4 

1 



x 



/ \ —a s / \ —a t 



(1 + S\)(l + s|)(l + S3) I S43 1 1 S32 1 1^21 1 |^4i I \ S 3lS42/ \ S 3lS42/ 



(E.l) 



where we have introduced the notation 



(E.2) 
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In Eq. (E.l) the integration is only over s\, s 2 and S3, while S4 is chosen in a way to 
preserve the cyclic symmetry but otherwise arbitrary. The difficulty in analyzing the 
amplitude (E.l) is the presence of the factors 1/(1 + s 2 ) which violate the PSL(2;M) 
symmetry. 

Let us introduce the variable 



x 



S31S42 

which runs from at s 2 = Si to 1 at s 2 = s 3 and vise versa. We can use x instead of S2: 

S41S43 



s 2 = s 4 - 



and 



ds 2 



■S43 + 

S43 s 31 s 41 dx 



(1 + si) s\ 3 (l + s\) + 2a;s43S3i(l + S1S4) + x 2 s 2 l3 (l + s 2 4 ) 
For the amplitude (E.l) we get 



(E.3) 

>f s 2 : 
(E.4) 

(E.5) 



G 4 



&xx- a ' s - l (l-xY alt - 1 



1 



1 



X 



S4 S3 S4 +00 +00 S3 

ds 3 J dsi + y ds 3 y dsi + y ds 3 y dsi 

— OO —OO —OO S4 S4 S4 

(s 43 + XS31) 2 



(1 + S 2 )(l + si) |s 43 ||s3i||s 4 i| [s| 3 (l + S 2 ) + 2xS 43 S 3 i(l + S1S4) +x 2 sf 3 (l + s|)] ' 

(E.6) 

where s 4 is not necessarily = 00, so the ordering of the points Si, s 3 and s 4 preserves 
cycling symmetry. 

The integrand in Eq. (E.6) differs from the projective-covariant one only by the ugly 
factor 



1 



(S43 + xs 31i 



(1 + s?)(l + s|) [4,(1 + sj) + 2xs 43 s 3l (l + Sl s 4 ) + x 2 s 2 3 (l + s 2 )] (K7) 
which we shall now see is not important because the integral over si and s 3 is divergent. 
In order to regularize the integral, we proceed like in Eq. (4.5) changing 



sj I + 1 - b , 



(E.8) 



I 2 J j I I °z °j I 

and rescale Sj — > (1 — The linearly divergent part of the integral over si and s 3 in 
Eq. (E.6) decouples from the integral over x and reads 

§4 S3 S4 +OO +OO S3 

- — - y ds 3 J ds! + y ds 3 J ds 1 + J ds 3 J ds x 



—00 —00 —00 S4 
1 



s 4 



s 4 



, ^ , ^ , ^ + 0((l-b) ). (E.9) 

(|5 43 | + l)(|53i| + l)(|54i| + l) lV >' K J 

This integral is convergent and does not depend on s 4 because of the invariance under 
translations. 

We have thus reproduced Eq. (4.16), modulo a constant given by Eq. (E.9), with the 
straight Regge trajectory a(t) = a't. The Regge behavior comes about in the same way 
as when there is the PSL(2; M) symmetry. 
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